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Abstract 

We incorporate realistic pair potential energies directly into a non-linear least-square fit of diffraction data to quanti- 
tatively compare structure models with experiment for the Ni-rich d(AlNiCo) quasicrystal. The initial structure models 
are derived from a few a priori assumptions (gross features of the Patterson function) and the pair potentials. In place 
of the common hyperspace approach to the structure refinement of quasicrystals, we use a real-space tile decoration 
scheme, which does not rely on strict quasiperiodicity, and makes it easy to enforce sensible local arrangements of 
the atoms. Inclusion of the energies provides information complementary to the diffraction data and protects the fit 
procedure from converging on spurious solutions. The method pinpoints sites which are likely to break the symmetry 
of their local environment. 
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1. Introduction 

The decagonal AINiCo system exhibits a variety of 
metastable and stable phases, out of which the "basic 
Ni-rich" phase has attracted many detailed stud- 
ies due to its high structural quality. Recently, two 
high-quality refinements have been published based 
on single-crystal X-ray E] ■ At the same time, a 
Monte Carlo method was developed 0] to predict the 
same structure, approximating the Hamiltonian with 
pair potentials, and using as inputs only the symme- 
try, the quasilattice constant and the fact that the 
structure is strongly layered. Here, we present the 
first trial of a combined method which uses energy 
and diffraction information to fix complementary fea- 
tures of the atomic structure. 

The refined structures of [2] and |3] agree on basic 
features, but they differ in many important details. 
The refinement presented by CHS [2] exhibits 8A pe- 
riodicity along the vertical axis, lacks a 10-fold screw 
symmetry property, and does not lend itself an easy 
interpretation in terms of a tiling geometry. The 
structure solution of |3j has 4A periodicity, and can 
be readily interpreted in terms of the Hexagon-Boat- 
Star (HBS) tiling with edge length of a q ~6.5A. The 



R-factors achieved are comparable in the two cases. 

Quasicrystals, as aperiodic structures, typically in- 
clude local patterns which are rare enough to be ill- 
determined from diffraction, but common enough to 
have a significant effect on physical properties (such 
as total energy or conductivity). Thus the problem 
of structure prediction from energies becomes entan- 
gled with structure fitting from diffraction. We sus- 
pect that, in the future, the best quasicrystal struc- 
ture fits will combine the two inputs in some fashion. 

Here, we present the first structure fit in quasicrys- 
tals, in which energy and diffraction information are 
combined in the same objective function. Previously, 
powder diffraction data on z-TiZnNi jS] was system- 
atically refined with energy inputs. In that work, 
each type of information was used separately to fix 
the degrees of freedom which are most sensitive to 
it: a least-squares refinement of diffraction data that 
determined the species occupying each site type was 
alternated with ab-initio calculations to relax the po- 
sitions. 

In the present study we calculate energetics from 
pair potentials microscopically derived using "gen- 
eralized pseudopotential theory" a total of six 
functions (e.g. VAiNi(r)) for all combinations of the 



species. The prominent feature of these potentials 
is Friedel oscillations, so that a typical potential 
has three minima before the cutoff radius (which we 
chose to be 10 A). The atom sites are called {R}, 
with x(R) designating the chemistry (Al, Ni, or Co) 
of that site. A convenient diagnostic tool in the re- 
finement of a decoration structure [7] is the "site en- 
ergy" , the portion of total energy ascribed to inter- 
actions of site R: 

^ itc (R)^^V; (g)x(fr) (|R-R'|). (1) 

R' 

Then the total energy is 

£ tot ^iX> site (R). (2) 

RR' 

Whereas previous refinements of quasicrystals 
have used the hyperspace-cut representation of the 
structure, our results demonstrate the equal ef- 
fectiveness of an alternative approach, the tile- 
decoration framework. The value of the tiling- 
decoration approach is that the structure models do 
not rely on periodicity or quasiperiodicity, and they 
are naturally represented in "real" space. Conse- 
quently, decoration is an elegant way [5] to trans- 
fer/combine information between the different sys- 
tem sizes appropriate for different calculations, e.g.: 
(i) numerical diffraction calculations from a large ap- 
proximant; (ii) eventual Monte Carlo simulations us- 
ing a tile reshuffling update and a tile-tile Hamilto- 
nian, in a smaller approximant; (iii) molecular dy- 
namics using pair potentials; (iv) atom relaxation 
using an ab-initio total energy code, which is typi- 
cally feasible only in the smallest approximants. 

The decoration formulation might be attractive for 
dealing with randomness (an issue we ignore in the 
present contribution, however). A realistic structure 
model is an ensemble of atomic configurations. When 
the contents of different sites are independent ran- 
dom variables, then a stochastic, mixed occupation 
of each site is a fair description. If, however, the 
sites are highly correlated, a random tiling may be 
a more reasonable model; indeed, that is our present 
picture of the "basic Ni" phase of d-AlNiCo 0j, A 
tile-decoration is superior to a hyperspace descrip- 
tion in the random-tiling case: all the randomness 
may be ascribed to the tile configurations, while the 
decoration of each tiling is deterministic. Some of 
the technical questions of a random-tiling structure 
refinement were discussed in Sec. 5 of Ref. [8]. 

2. Structure modeling for combined fit 

To explain the decoration approach, we review 
nomenclature introduced in [jj: A decoration is a 
mapping which, given any valid tile-configuration, 



produces a set of atom sites. The sites are grouped 
into orbits, each of which is bound to a particular 
kind of tile, and lies in the same positions on every 
tile of that kind. An orbit plays a role like an orbit 
of symmetry-related Wyckoff positions in a crystal 
structure. In the hyperspace-cut representation of a 
quasiperiodic structure, each member of the decora- 
tion orbit on a representative tile would correspond 
to one subdomain of the acceptance domain in per- 
pendicular space. 

More generally, orbits are bound to tiling objects 
which may be not only tiles, but other geometrical 
constituents of the tiling such as tile vertices - e.g. 
we may choose to surround every vertex by an identi- 
cal symmetric cluster. Another kind of a tiling object 
is defined to overlay a certain local pattern in the 
tiling. We can then introduce context dependence 
by rebinding certain atoms to the new tiling object. 
Other atoms remain bound to the old objects, so 
we do not define too many positional parameters. 
One decoration rule is a rebinding of another if it 
has more orbits, and the family of structures pro- 
duced from the original decoration (by varying the 
parameter values) is a subset of the family produced 
from the rebinding decoration. 

By default we require that atomic decorations in- 
herit the symmetry of the tiling objects so that plac- 
ing an atom at general position will generate all sym- 
metry equivalent sites. Furthermore, preservation of 
the site symmetry constrains the degrees of freedom 
for atomic displacements. In reality, forcing atomic 
motifs to obey the point symmetry of tiling object is 
not always energetically favourable. 

In the tile-decoration machinery, such situations 
are handled by the concept of symmetry-breaking: a 
flavour (like an arrow on a tile) is associated with 
each (broken) symmetry related variant. |7J, so that 
the decoration remains formally deterministic and 
reflects the symmetries of the flavoured tiles. The 
flavours may be initialized randomly, and in the sim- 
ulation process (for example Monte Carlo annealing) 
the symmetry-breaking may be effectively used as an 
update move. 

The decoration rule (see FigQJ for our combined 
fit is adapted from the deterministic model of U, in 
which atoms decorate an HBS tiling with edge length 
a q ~6.4A(similar to Fig. 3a of 0], but without the 
complication of the "bow-tie" defects shown there.) 
Atoms lie on flat layers that are spaced by 2A and 
are related by the IO5 screw axis. Pairs of adjacent 
flat layers constitute a "4A bilayer" which can be 
taken as the fundamental repeat unit of a 4A periodic 
structure but can also be stacked in other fashions. 
This model was derived by successive elimination of 
the degrees of freedom in Monte-Carlo simulations 
with the GPT pair potentials. 

Notice that our Hexagon tile has a mirror plane 
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running down its long axis, but not transverse be- 
cause the Al(9) atom breaks this symmetry. We 
should visualize it as if marked by an arrow to dis- 
tinguish the two ends. In the same way, the Boat's 
mirror symmetry and the Star tile's fivefold rota- 
tion symmetry are broken (respectively by Al(12) 
and Al(16) atoms) to allow for our decoration. 

Our "basic" rule uses as supplementary tiling ob- 
jects the Fat rhombi (which appear naturally when 
the HBS tiles are subdivided into Penrose rhombi), 
as well as the interior vertices thereby created; this 
captures the similarity of the local environments in 
all the 72° corners of HBS tiles. This rule is shown in 
Figffl All Ni atoms are ascribed to the HBS edges, 
Co atoms decorate the supplementary interior ver- 
tices, as well as a site on the long diagonal of every 
Fat rhombus. Orbits 12 and 16 implement the bro- 
ken symmetries of the Boat and Star. In the end, this 
decoration rule associates 16 orbits with nine tiling 
objects, and yields 38 positional parameters. 

Refinement using the "basic" binding motivated a 
more complex "augmented" binding (Table [JJ, that 
proved to be more successful in fitting the data. 
Firstly, since the Ni(3) orbit alone contains more 
than 20% of all atoms, we split it among 3 objects, 
"context-sensitive" HBS Nodes. In both approxi- 
mants we used, three classes of HBS vertex stars oc- 
cure: 5-fold, 3-fold and 2-fold (in random HBS, also 
4- fold and another 3- fold vertex type occurs). The 
symmetry of these objects gives rise to 4 Ni orbits, 
labelled Ni(3 a -3c;). The new orbits are listed in Table 
^ below the horizontal line. Second, one of the two 
Al(15) atoms has been replaced by a vacancy, imple- 
menting symmetry-breaking of the Hexagon mirror 
(Al(15+) in TableUJ), and new atom Al(17) has been 
placed near the center of the Hexagon. 

The atomic density for the decoration rule (in ei- 
ther binding variant) is 0.0705 atoms/A 3 , and the 
composition is Al70Ni20.7C09.3j slightly richer in Co 
than the CHS sample (Al70.6Ni22.7C06. 7). When 
the atoms are fixed at "ideal" sites, the minimal 
distance occuring between Al-Al, Al-Co and Al-Ni 
pairs is 2.46A. The structural energy of this initial 
model is —0.297 eV/atom for a small approximant 
with 214 atoms per unit cell; using positional de- 
grees of freedom of the decoration rule it relaxes to 
—0.393 eV/atom. Fully relaxed energy (atomic dis- 
placements for atoms in the same orbit no longer 
equivalent) is —0.422 eV/atom, and finally allowing 
atoms to relax their positions in a double period (8A) 
unit cell, the energy drops to —0.442 e V/atom. 

A drawback of the tile-decoration approach is that 
the infinite non-periodic quasicrystal must usually be 
handled indirectly, being approximated by an peri- 
odic "approximant" . This necessarily has a "back- 
ground phason strain" which is, roughly, the average 
tilt of the cut through hyperspace away from the 



irrational orientation with 10-fold symmetry.) Con- 
sequently, each orbit of reciprocal lattice vectors Q 
of the quasicrystal gets split into several incquivalcnt 
orbits of the approximant 's reciprocal lattice vectors 
Qapp- In the fit, we selected only one representative 
Qapp for each decagonal Q, since the number of dat- 
apoints is already large, but we monitored system- 
atically the errors thus introduced, calculating <7 app 
as the r.m.s. deviation of the difraction amplitudes 
among all Q app mapping on a given Q. For the larger 
approximant discussed below, <r app was similar to or 
less than the experimentally observed uncertainty for 
all peaks. 

A more systematic treatment of the approximant 
diffraction amplitudes would be to symmetrize over 
all those deriving from a single orbit of quasicrystal 
Q. This would demand careful handling of the phase 
factors, by properly centering the approximant struc- 
ture both in real and in perpendicular space 0- 

A further difficulty is that for smaller approxi- 
mants, different quasiperiodic Q vectors may map on 
the same approximant Q app . We found the smallest 
orthorhombic approximant avoiding such ambiguous 
mappings for all reflections in the CHS data set has 
edge lengths 6lA x 32A in the decagonal plane. In 
the following, we label this approximant T42 (for the 
number of HBS vertices per unit cell). For our dec- 
oration rule it contains 560 atoms per unit cell. We 
found it satisfactory to use just one more bigger ap- 
proximant, with both sides in the decagonal plane 
larger by a factor of r = (y/E+ l)/2. We denote this 
tiling Tno and decorate it with 1466 atoms. 

3. Procedure 

We used the CHS dataset • which contains 2767 
unique reflections under the assumption of 10/m 
Laue group; by averaging over symmetry-related am- 
plitudes, to enforce 10/mmm Laue symmetry, a set 
with 1544 unique reflections is obtained. The two 
possibilities gave about equal internal i?-factors. Our 
structure model has PlO^/mmc space group, so we 
used the smaller 10/mmm data set. (Using the 10/m 
set gave practically the same results, except that typ- 
ically both the K w and R factor are increased by 1%.) 

The partial diffraction amplitude for tiling orbit 
j is 

F° lh (Q) = J2w ja (Q)UQ) E e ^'* ( 3 ) 

The inner summation in Q runs over all atoms of 
species a in orbit j, and f a (Q) is the atomic form 
factor for that species. The Debye- Waller (DW) coef- 
ficients b xy j a and b z j <a couple with decagonal-plane 
and vertical Q components, and take different values 
not only for different orbits, but also for different 
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species occupying the same orbit: 

W ja @) = exp[~(b*y ja §l y + b z ja Ql] (4) 
The total calculated diffraction amplitude is then 

^ ca/c (Q) = e-4 6j -3" \F° Th \j{Q)- (5) 



where we introduced a single perp-space DW coef- 
ficient b±. Finally, we compute goodness/reliability 
factors 



,2 _ 



X 



^(F ofcs (Q)-F ca/c (Q)) 2 /a(Q) 2 (6) 



F obs (Q) - F ca ' c (Q)| 



(7) 



As usual, the weighted i?-factor, called R w , is defined 
like (JJJ except that each term in the numerator and 
denominator acquires a factor l/er(Q) 2 . 

Implementation. We couple pair-potential energy 
with x 2 v i a parameter A in the objective function for 
the non-linear least-squares minimization: 



U = x + [E - E tars Y/X 



(8) 



where E is the energy J3J), per atom. Here E tai ' s is 
a "target" energy, 1 and is treated in the fit as a 
dummy datapoint with an "error bar" A. E tars may 
be set to the total structural energy per atom after 
an (unconstrained) relaxation E™t. 

Our fitting form (Eq. allows independent DW 
parameters b xy and b z for each orbit/species combi- 
nation, but these are underdetermined by the data; 
to obtain sensible results, the fit needs to be biased 
towards having similar DW factors for similar atoms, 
in the spirit of the maximum-entropy method. We 
implement this by adding terms for b xy and b z to 
the x 2 sum: before each iteration, we calculate the 
average DW factor for each species and the xy and 
z components of DW factor, and set them as data- 
points for the next iteration, with appropriately cho- 
sen <tdw- 

4. Results 

Fitting the "basic" binding to the CHS dataset we 
obtained reasonable, but not satisfactory R-factors 
R w ={).\22 and #=0.202. The large displacement 
of atom Al(15) (see Fig. El top) and other con- 
siderations motivated our "augmented" binding; Ta- 
ble El summarizes our results for the latter. The 

iWhen AE = E tot - £ tar s > A, as is usually the case in 
the orbit-constrained relaxation, then U = x 2 +2(A_E/A 2 )E tot 
would be practically equivalent to JSJ; since the likelihood of 
parameters is supposed to behave as exp(— U), X 2 /2AE is 
evidently playing the role of a temperature. 



symmetry-constrained relaxed energy E Te \ of the 
"augmented" rule is tolerably higher than E IC \ of the 
"basic" rule (see TableEJ). However, the combined fit 
improved dramatically, displacing the Al(15 a ) atom 
towards the long body diagonal of H tile. We have 
selected two structure solutions by the combined fit, 
one with lowest R^, factor but rather high energy, 
and one with slightly worse R^, factor, but compara- 
ble in energy with the basic decoration rule (Table 

EJi. 

The F ca i c (Q) from our fit is compared to the 
F bs(Q) from the diffraction data in Fig. EJ 

The fitted DW factors, especially those of Al, ex- 
hibit high degree of anisotropy (Table El )■ This 
agrees with the TYT refinement [3] of a different 
sample of AINiCo. (The TYT sample was slightly 
richer in Al and Co than the sample for the CHS 
data that we are fitting here.) 

Most of the TYT hyperatom orbits correspond 
to our HBS tile decoration orbits, and their indi- 
vidual orbit DW factors are mostly consistent with 
ours for the corresponding atoms. The most pro- 
nounced disagreement is just the orbit Al(15+) for 
which b z /b xy « 5 in our refinement, but sa 1/5 in 
the TYT refinement (see Table 2 in Ref. |3J, site 
label 9). This atom is at the same time displaced 
farthest from its ideal position for both TYT and 
our refinements, and the displacements have simi- 
lar magnitude and direction for TYT and our "ba- 
sic" binding (see FigEJ). In the "augmented" bind- 
ing, Al(15+) occurs in a completely different envi- 
ronment, displacing towards long diagonal of H tile, 
which explains the discrepancy in the DW factors. 

The TYT model and our decoration (in either vari- 
ant) differ from the CHS model in having a larger 
space group symmetry (10-fold screw axis vs. a 5- 
fold axis, 4A translation vs. 8A translation), but also 
in other respects: higher atomic density, and absence 
of transition metals (TM) in the sites correspond- 
ing to hyperatom "B" in Ref. Our decorations 
are more constrained than the other models, in that 
mixed Al-TM occupancy is disallowed, partial oc- 
cupation appears only through symmetry-breaking, 
and unphysical close pairs never appear on account of 
the energy term in Eq.JSJl The CHS model achieved 
a lower weighted R-factor (^0.06) than the present 
refinement of the same data set, at the price of in- 
cluding many more parameters. Enlarging the set 
of fitting parameters within our combined energy- 
diffraction approach should further reduce the R- 
factor while maintaining our energy-based tile deco- 
ration description. 

5. Discussion 

The rearrangement of the atomic decoration inside 
the Hexagon tile - driven, we think, by the diffrac- 
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tion term in (jSJ) - is reminiscent of the atomic shifts 
observed in molecular dynamics (MD) simulations 
under the same potentials with an initial atomic 
configuration similar to the present model Kj|. In 
those simulations, using supercells 8A or 16A thick 
in the stacking direction, slow cooling led to a more 
favourable arrangement of atoms, in which one Al 
was pushed from the Al(15) position towards the cen- 
ter of the Hexagon, while the other moved onto the 
long body diagonal of the Hexagon. However, in con- 
trast to our present model in which two such atoms 
get squeezed into a column, in 9 only one atom 
per 8A period was displaced in this fashion, so that 
the column in the Hexagon center contained three 
Al atoms per 8A period. A recent MD study [TO] 
confirms that the optimal Hexagon decoration has 
8A period, with atomic arrangement topologically 
equivalent to the result of our combined fit, but in 
[TU| the atomic sites in the second 4A bilayer are 
flipped in each Hexagon across a mirror plane per- 
pendicular to its body diagonal. 

These findings again call for study of the so-called 
"stacking-fault" mechanism, in which the decoration 
rule would remain essentially a 4A rule, and the vari- 
ations in the decoration from one 4A bilayer to the 
next would be described by "flips" of the tiles, oc- 
curing independently in each bilayer. 

The discussed variations in the atomic structure 
also illuminate a motif generally considered impor- 
tant in decagonal structures: the 20A diameter 
decagon. One such cluster in marked in Fig. by 
the inscribed circle. Note the inner ring containing 
~20 Al atoms (Al(17) and Al(9)) "squeezed" in be- 
tween 10 TM atoms (typically Co, sometimes Ni), a 
feature that emerges in many HREM/HAADF im- 
ages ^]^]. This is energetically unfavorable in a 
structure with strict 4 A periodicity due to crowding 
of the repulsive Al atoms. It occurs in our combined 
fit for a 4A bilayer only because it is favored by the 
diffraction data. On the other hand, as suggested 
by the MD studies, it should emerge as a low energy 
configuration from simulations performed for stacked 
bilayers which relax the strict 4 A periodicity. 

In conclusion, we used tile-decoration machinery 
to set up a combined energy-diffraction data fit of 
the d-AlNiCo structure. The resulting R-factors are 
comparable to those from a previous study 0, while 
using fewer parameters, and not allowing ad-hoc av- 
eraged/mixed site occupancies. Apart from the in- 
clusion of energies into the diffraction fit, at this 
stage our refinement is completely analogous to the 
hyperspace refinements, since we used strictly ap- 
proximants of the ideal, unperturbed quasiperiodic 
tiling. The possible disorder correlations hidden in 
the experimental data are still lumped into global pa- 
rameters of the fit, mainly the perp-space DW factor. 

A future prospect is to take full advantage of 



the real-space formulation of the problem and re- 
cover the correlations in atomic occupancies directly 
from a "grand combined fit" . This would use Monte 
Carlo annealing to sample all relevant degrees of free- 
dom (swapping individual atoms, flipping symmetry- 
broken tiling objects) in conjunction with optimizing 
R-factors fitting the diffraction data. |Sj If the pair 
potentials prove indeed sufficiently realistic, the fit 
parameters may be further constrained by DW fac- 
tors pre-calculated from zero-T phonon spectra |13j . 
or extracted from room temperature MD annealing 
runs. 

Finally, we expect our approach will prove opti- 
mal for refinements of large quasicrystal approximant 
structures, - several such "approximant" phases are 
known to exist in the AINiCo system - that 
are too complex to be handled by routine crystal- 
lographic refinement approaches. 
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Table 1. Atomic orbits in "augmented" binding, x is 
fractional content, b xy and b z refined anisotropic DW co- 
efficients, and E S it e are site energies in eV/atom, for the 
approximant Tno- Orbits above the horizontal line were 
unchanged from the "basic" binding. Superscript r de- 
notes objects defined from the corners of the vertices in- 
side the HBS tiles (Fat r and Bond r ), see also Fig. 
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Fig. 1. Atomic decoration of the HBS tiling for "ba- 
sic" (top panel) and "augmented" (bottom panel) bind- 
ings. Dark circles are Ni atoms, dark triangles Co, gray 
Al. Unoccupied sites are empty circles. Tiling orbits 
are shown as framed labels (once for each tiling object). 
Thick lines are HBS-tiling edges, thin lines edges of the 
inscribed rhombus tiling. For orbits that are binded to 
HBS tiles, the edges of rhombuses inside the tiles are not 
drawn. Each symbol occurs in two sizes, for upper and 
lower layer in the 4A thick slice cut perpendicularly to 
the periodic direction. 
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binding 
basic 
augmented 



9 
13 



iV orb 


Npos 


iV par 




16 


38 


72 


-0.400 


23 


58 


106 


-0.395 


bperp 


R w 


R 


E 


3.65 


0.122 


0.202 


-0.372 


4.15 


0.089 


0.169 


-0.344 


3.88 


0.082 


0.152 


-0.330 


3.93 


0.086 


0.159 


-0.370 



model 
T4 2 -basic 
T42-augmented 
Tuo-augmented 
(smaller A) 



Table 2. Top: Numbers of tiling objects iV bj, orbits 
A^orb, positional parameters N pa s and total number of 
fit parameters N paT (including anisotropic DW factors 
and the global b± DW factor). The last column is the 
energy per atom after constrained relaxation, without use 
of diffraction data. Bottom: summary of the results for 
the combined fits, giving the perp-space Debye-Waller 
coefficient, the weighted (Rw) and unweighted R factors. 
The last row reports an alternative refinement result for 
approximant Tno with the "augmented" decoration rule, 
in which the A parameter in Eq. |S]was kept smaller to 
achieve lower energy. 




Fig. 3. Observed vs calculated diffraction amplitudes 
for approximant Tno, final refined model of the "aug- 
mented" binding. 



Fig. 2. Top: "basic" decoration rule, refined model. 
Arrows indicate displacements (empty symbols stay at 
the "ideal" starting position) that occured during com- 
bined fit. Vacant sites are shown as small circles and 
marked by crosses. Bottom: "augmented" binding, re- 
fined model. One decagonal columnar cluster is marked 
by two dashed circles, near the outer shell of atoms, and 
near the inner ring of 10 TM atoms. 
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